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1. Introduction 

The main analytical method used to treat the physics of interacting systems is 
perturbation theory. In this framework, bare, noninteracting, matter fields and gauge 
fields (to which bare quanta are associated) are introduced, the interaction being 
accounted for order by order in powers of the coupling constant in the perturbative 
expansion of the observables. For a long time, physicists have been aware that 
perturbative expansions have a limited scope, in spite of all the remarkable achievements 
obtained with them. Some examples of such situations are resonant effects associated 
to the coupling of atoms to strong radiofrequency fields, or the low-energy domain 
of quantum chromodynamics, where confinement of quarks and gluons takes place. 
Along the last decades, many attempts were devised to circumvent the limitations of 
perturbation theory, in particular when strong effective couplings are involved. 

In any case, as a matter of principle, due to the nonvanishing of the coupling 
constant, the idea of a bare particle associated to a bare matter field is actually an 
artifact of perturbation theory and strictly speaking is physically meaningless. In 
general, a charged physical particle is always coupled to a gauge field; in other words, 
it is always "dressed" by a cloud of quanta of the gauge field, for instance, photons in 
the case of electrodynamics. 

With respect to the time evolution of systems of (matter) particles, the idea is that 
particles are coupled to an environment or to a thermal bath. There are usually two 
equivalent ways of modeling the environment to which the particle is coupled: either to 
represent it by a free field, as was done in Refs. [H [2], or to consider the environment 
as a reservoir composed of a large number of noninteracting harmonic oscillators (see, 
for instance, P HI El E] ) . In both cases, exactly the same type of argument given above 
for a charged particle applies, making the appropriate changes, to these systems. We 
may speak of a "dressing" of the set of particles by the ensemble of the harmonic modes 
of the environment. The particles in the system are considered as always "dressed" by 
a cloud of quanta of the environment. This is true in general for any system in which 
material particles are coupled to an environment, no matter what the specific nature of 
the environment and interaction involved are. 

In the present work we adopt this point of view and consider a harmonic oscillator, 
coupled linearly to an environment modeled by the infinite set of harmonic modes 
of a scalar field, the whole system being contained in a perfectly reflecting sphere 
of radius R. We implement a nonperturbative study of the time evolution of a 
superposition of the first excited and the ground states of the oscillator, by means of the 
dressed states introduced in [7] and already employed to investigate several situations 
[HI [101 [HI [121 US] • Such a nonperturbative approach is possible due to the linear 
character of the interaction between the oscillator and the environment we choose. 

The semiqualitative idea of a "dressed atom" was originally introduced in Refs. [H] 
and [15], and largely employed in studies involving the interaction of atoms and 
electromagnetic fields [T8l [T9l [201 l2l] . In the realm of general physics, the concept 
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of dressing a matter particle by an environment has found an application in describing 
the radiation damping of classical systems [16]. Our dressed states can be viewed as a 
rigorous version of these dressing procedures. Moreover, the oscillator may represent a 
mode of the quantized electromagnetic field in a cavity interacting with the environment; 
in this case, our study refers to the time evolution of the superposition of the ground 
state and the first excited state of the field in the cavity. 

The model in itself is not new; a system composed of a material body interacting 
linearly with an environment has been the main subject of many papers, as those quoted 
in Refs. [H [21 El IH O El [T7l [16] . The novelty lies in the nonperturbative approach to the 
problem, by means of renormalized coordinates and dressed states, which was started in 
[7]. In particular, our dressed states are not the same as those employed in the literature, 
usually associated to normal coordinates. Our dressed states are given in terms of our 
renormalized coordinates and allow a rigorous study of the time evolution of quantum 
systems. The results we obtain by this means are those expected on physical grounds, 
but contain corrections with respect to the formulas obtained from perturbation theory. 

We structure the paper as follows. In Section II we review the basic aspects of the 
formalism of renormalized coordinates and dressed states and then, in Section III, we 
treat the time evolution of a superposition of oscillator states in both limiting situations 
of R — > oo and small R. Section IV contains our final comments. 



2. Renormalized coordinates and dressed states 



Let us start by considering the harmonic oscillator having bare frequency ujq, linearly 
coupled to a field described by N (— > oo) other oscillators, with frequencies u>k, 
k = 1, 2, . . . , N. The whole system is contained in a perfectly reflecting spherical cavity 
of radius R, the free space corresponding to the limit R —>■ oo. Hereafter, we shall refer 
to the harmonic oscillator as the particle, to distinguish it from the harmonic modes of 
the environment. Denoting by qo(t) {po(t)) and quif) ij>kif)) the coordinates (momenta) 
associated with the particle and the field oscillators respectively, the Hamiltonian of the 
system is 



H 



N 



PC + ^0<?0 + E ( P k + U k<ll) 



k=l 



N 



qo^2r]UJ k q k , 



(1) 



k=l 



where the limit N — > oo is understood and rj is a constant. 



The Hamiltonian 
transformation, 



JT]) can be turned to principal axis by means of a point 



N 



A' 



Ew 



E£ p - 



(2) 



r=0 



where \i = (0, {&}), k = 1,2, 
The subscripts \i = and fi - 



.,N, performed by an orthonormal matrix T = (t r ) 
k refer respectively to the particle and the harmonic 



Time evolution of a superposition of dressed oscillator states in a cavity 



4 



modes of the field and r refers to the normal modes. In terms of normal momenta and 
coordinates, the transformed Hamiltonian reads 



1 N 



2 ) 



(3) 



where the Q r 's are the normal frequencies corresponding to the collective stable 
oscillation modes of the coupled system. 



Using the coordinate transformation ([ 
making use of the normalization condition 

JV 

£(Q 2 = i, 

Ai=0 



in the equations of motion and explicitly 



(4) 



we get 



f 



JV 



v 2 ^l 



fe=i 



with the condition 



E 



fc =i ^ 



nr 



The right hand side of Eq. (jSJ) diverges in the limit N 
5u 2 = Nrj 2 , it can be rewritten in the form 



N 1 

5u 2 - n 2 r = v ^ r -r 1 : 



fc=l k 



n 2 ' 



(5) 



(6) 

oo. Defining the counterterm 



(7) 



Equation ([7D has + 1 solutions, corresponding to the iV + 1 normal collective 
modes. It can be shown [7] that if u 2 > 5oo 2 , all possible solutions for Q 2 are positive, 
physically meaning that the system oscillates harmonically in all its modes. On the other 
hand, when u 2 < 5u 2 , one of the solutions is negative and so no stationary configuration 
is allowed. Nevertheless, in a different context, it is precisely this runaway solution that 
is related to the existence of a bound state in the Lee-Friedrichs model. This solution 
is considered in Ref. [22] in the framework of a model to describe qualitatively the 
existence of bound states in particle physics. 

Therefore, we just consider the situation in which all normal modes are harmonic, 
which corresponds to the first case above, Uq > Sou 2 , and define the renormalized 
frequency 



u 2 



lim (uq - Nrj 2 ) 



N^oo 



following the pioneering work of Ref. [23] . In the limit iV 

V ^u 2 -W 



oo, Eq. (|7j) becomes 



CO 



(8) 



(9) 



k=l 



We see that, in this limit, the above procedure is exactly the analogue of mass 
renormalization in quantum field theory: the addition of a counterterm —Nrfq 2 allows 
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one to compensate the infinity of Uq in such a way as to leave a finite, physically 
meaningful, renormalized frequency to. 

To proceed, we take the constant r\ as 



V 



4gAu 



71 



(10) 



where Au; is the interval between two neighboring field frequencies and g is the coupling 
constant with dimension of frequency. The environment frequencies u k can be written 
in the form 

7TC 



u k = k- 



R 



k = 1,2, 



and, so, Au = txc/R. Then, using the identity 



E 

k=l 



k 2 



u 



1 

2 ~ 2 



7T 



u 



cot (ttu) 



Eq. can be written in closed form: 

ft c 

cot {— )=Y g + m 



-,2 



2^ 



(12) 



(13) 



The elements of the transformation matrix, turning the particle-field system to principal 
axis, are obtained in terms of the physically meaningful quantities ft r , u, after some 
rather long but straightforward manipulations [7J. They read 



(ft 2 - Q 2 ) 2 + \ (3ft 2 , - Q 2 ) + 4g 2 n 2 



,r _ -,-.» 



(14) 



(15) 



Let us now consider the eigenstates of our system, |n , rii, n 2 , ...), represented by 
the normalized eigenfunctions, written in terms of the normal coordinates {Q r }, 



ih..xQ>t) = n 




ft. 




r n e 



(16) 



where Hi s stands for the l s -th Hermite polynomial and T is the normalized vacuum 
eigenfunction. 

We introduce renormalized coordinates q' and {q^} for the dressed particle and 
the dressed field, respectively, defined by 

V^ = X)*MV^n (17) 

r 

where a) M = {u, Ui}. In terms of the renormalized coordinates, we define for a fixed 
instant, t = 0, dressed states, \kq, ki, k 2 , ...) by means of the complete orthonormal set 
of functions [7] 



'2 k m 
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where q'^ = {q' Q , g-}. Note that the ground state T in the above equation is the same as 
in Eq. (fT6l) . The invariance of the ground state is due to our definition of renormalized 
coordinates given by Eq. ( IT7I) . Each function ip KoKl ...(q') describes a state in which the 
dressed oscillator q'^ is in its K M -th excited state. In terms of the bare coordinates, the 
renormalized coordinates are expressed as 

^ = ^2 a ^9^ ( 19 ) 

where 



E^V^- (20) 



Remark that the introduction of the renormalized coordinates implies, differently 
from the bare vacuum, the stability of the dressed vacuum state since, by construction, 
it is identical to the ground state of the interacting Hamiltonian ([3]). Also it is important 
to notice that the renormalized coordinates and the dressed states, defined in Eqs. (fT7j) 
and (118p are new collective objects, different from the normal coordinates Q and the 
eigenstates ffl6|) . Since the transformation (fTTj) is not orthogonal, the Hamiltonian is 
not diagonal in the renormalized coordinates. Thus, distinctly from the eigenstates, our 
dressed states are all unstable, except for the ground dressed = 0}) state. We shall 
assume that the dressed states (fT8l) are the physically appropriate states to describe 
the time evolution of superpositions of states of the mechanical oscillator, taking into 
account non-perturbatively the effect of the interaction with the environment. This is 
an alternative to the use of states written in terms of the bare coordinates g M , which 
would require a perturbative renormalization procedure to correct order by order the 
oscillator frequency. 

Let us consider, for instance, the particular dressed state |r^(0)), represented by 
the wave function ^00-1(^)0- (?')■ ^ can be seen as describing, at a given instant, the 
configuration in which only the /i-th dressed oscillator is in the "first excited level" , all 
others being in their "ground states". These "levels" should not be confused with the 
stationary states given by Eq. ( fl6l) ; from now on, we shall use such a terminology to 
facilitate the reference to dressed states. As shown in Ref . [7] , the time evolution of the 
state |r^) is given by 

IW)> = E^Wirr(o)>, (21) 



where 



Moreover, it can be shown that, for all /i, 

£lU*)| 2 = l, (23) 

V 

which allows to interpret the coefficients f^vit) as probability amplitudes; for example, 
foo(t) is the probability amplitude that, if the dressed particle is in the first excited state 
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at t = 0, it remains excited at time t, while foi{t) represents the probability amplitude 
that the i-th dressed harmonic mode of the cavity be at the first excited level. Also, 
the elements of the matrix TjP^'.., connecting the dressed states to the eigenstates, can 
be evaluated. For the dressed state in which only the K^-th dressed oscillator is in the 
N-th excited level, all other being in the ground state, we have [7j, 

^--(w^)*^ (24) 

where l + l\ + ■ ■ ■ = N. 

3. Behavior of the superposition of states in a cavity 

We now analyze the situation in which, at time t = 0, the environment oscillators are 
in their ground states and the particle is in a superposition of dressed states, and the 
whole system is allowed to evolve in time. Specifically, we take, as the initial particle 
dressed state, an arbitrary superposition of the first-excited and the ground states, 

|V>> = v^|r?(0)>+v^e*|r ), (25) 

where < £ < 1. While the ground state is stable, the particle first-excited 
state evolves in time according to Eq. ( |2T|) . that is, |r°(t)) = J2u fou(t) |r^(0)) with 
fov(t) = Yli S toK e ~ inat ■ At the instant t, the state of the system is given by the density 
matrix 

g(t) = e- iHt \^)(ij\e iHt . (26) 

We are interested in studying the influence of the environment on the time evolution 
of the state of the subsystem corresponding to the particle. Thus, we take the trace 
over all the degrees of freedom associated with the field, that is, we consider the reduced 
density matrix 

oo 

p(t)= (h,k 2 ,...\e- im \i;me im \k 1 ,k 2 ,...), (27) 

{k t =0} 

whose elements between dressed particle states are 

oo 

Pmn(t) = (m\p\n) = (m, h, k 2 , ...\e- iHt \ij)(tP\e iHt \n, k u k 2 , ...). (28) 

{k t =0} 

Replacing Eq. fTSBD in Eq. fT2"8"|) we get 

oo 

Pmn(t)= E {^ k2 --it)<^it)+ (l-O^00: fc2 --'W^000 lfc2 --*W 
{fei=0} 

+ \fW^V^^\t)<t '■"*(*)} , (29) 

where 

AZ°Z?. 2 r(t) = (m ,m 1 ,m 2 ,...\e- im \n ,n 1 ,n 2 ,...) (30) 
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are the probability amplitudes between the dressed states. These probability amplitudes 
can be computed using Eqs. fl22l) and fl24l) . We then have 

(31) 

A ^m 2 ... (t) = e -iE t /o/ ,(t)5 lmM J] 6om v , (32) 

where E = J2 r ^r/2 is the ground state energy of the system. Substituting the above 
equations into Eq. (1291) . we obtain, after summing over the fcj's, 



Pmnit) = - £ + £ ^ \fok(t)\ 2 j 5 0m 5 0n + £\foo{t)\ 2 5 lm 5 ln 

+ y/£(l-Qe- i +f 00 (t)5 lm 6 0n + ^(l-O^WWi, (33) 

(34) 

Using that ^ |/o/J 2 — 1 ; the nonzero elements of the reduced density matrix are given 
by 

Poo(i) = l-£|/ooW| 2 , (35) 
Pii(t)=t\foo(t)\ 2 , (36) 
PioW = y ^l-Q e-^/ooCt), (37) 
Poi(t) = V(M)^/;(^ (38) 

Note that Tr[p(t)] = 1, which means that the reduced density matrix p(t) does represent 
a physical state of the particle; also, it is not a pure state, since Tr[p 2 (t)] ^ 1. The degree 
of impurity of a quantum state can be quantified by the departure of the idempotent 
property; in the present case, we have 

D(t;0 = 1 -Tr[p 2 ] = 2e 2 |/oo(t)| 2 (l - |/oo(t)| 2 ) ■ (39) 

We are thus left with the calculation of the quantity foo{t), which will be done for the 
particular situations of an arbitrarily large cavity (R — > oo) and of a small cavity. 

3.1. Large cavity 

For an arbitrarily large value of R, Eq. (fT4l) reduces to 

f -> lim ^T^^TR m 

Also, in this limit, Au; = irc/R — > du = dQ and the expression for the quantity foo{t), 
given by Eq. (|22|) . can be cast in the form 

a poo Q 2 e~ int 

'»<«> = t/ " V-^ + vn* - («) 

The real part of foo(t) can be evaluated directly using the residue theorem. Defining 
the parameter k by 

k = v/^^ 2 , (42) 
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we obtain, for the three distinct cases, 
a) k 2 > : 



/oo(*) 



-gt 



cos nt — — sin nt 

K 



+ iG(t;u,g); 



b) k 2 = : 



f 00 (t) = e-*[l-gt]+iG(t ] u,g); 
and c) k 2 < : 

/oo(*) = 

where the function G(t; u, g) is given by 



-gt 



cosh lult — -, — r sinh IkH 

\k\ 



iG(t;u,g), 



G{t;uj,g) 



4g 



71 



dy 



y 2 sin yt 



(y 2 — uo 2 ) 2 + 4g 2 y 2 



(43) 
(44) 

(45) 
(46) 



The overall behavior of the function G(t; u, g) is illustrated in Fig. [TJ we see then that, in 
all cases, the real and imaginary parts of foo{t) decay with the time, for large times. This 
aspect dictates the behavior of the degree of purity, as a function of time, as illustrated 
in Fig. [21 for a case with g < u; similar results are obtained for the other regimes. 




Figure 1. Behavior of G(t; Q, g) as a function of t for the three distinct cases: Q = 1.5 
and g — 1.0 (full line); Q — 1.0 and g = 1.2 (dashed line); and uj = g = 2.0 (dotted 
line), in arbitrary units. 

Replacing / o(*), given by (gU) or (J45]), into Eqs. (I35jl - (j35]l leads to the 

elements of the reduced density matrix; simpler expressions are obtained for large times 
(t ^> 1/uj) at both weak- and strong- coupling regimes. For large t, G(t;u,g) can be 
approximated as 

Gft0,,)»J^ (47) 

Thus, in the limit case of weak coupling between the particle and the environment, 
g ^ uj (corresponding to k ~ u>), the large-time approximation of the reduced density 
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t 



1 2 3 4 5 6 



Figure 2. Behavior of D(t;£) as a function of t, with Q — 1.0 and g — 0.5 fixed (in 
arbitrary units), for some values of £: 0.3, 0.6 and 0.9 (dotted, dashed and full lines, 
respectively). 



matrix gives 



Pn(t) 
PooO) 
Pw{t) 
Poi(t) 



1 2 



l-Pll(t) 



COSCo>£ 



sin ujt 



-i<f> 



-<)t 



+ 



64g 2 1 

^2^6 J 



cos tut 



sin cut 



+ z 



8g ] 
nuH 3 J 



(48) 
(49) 
(50) 
(51) 



Similar results can be obtained for the case of a strong coupling between the particle 



and the environment, that is, when g ^> uj (i.e. \k\ 

64g 2 



g); for large t, we find 



Pu(t) 

Pw{t) 
Poi(t) 



e~ 4at + 



+ i 



8g 



TTUOH 3 



(52) 
(53) 
(54) 
(55) 



We see from the above equations that, for both weak- and strong-coupling regimes, 
in the case of a very large cavity, except for the vacuum-vacuum component which goes 
to 1, all other matrix elements of the reduced density matrix tend to zero as t —>■ oo. 
This means that the particle mixed state pit) tends to the dressed-particle vacuum state 
as time evolves; thus, for large cavities, the particle system is dissipative. This is also 
illustrated in Fig. El with £)(£;£) vanishing as t — > oo for all values of £. Notice that the 
the real and imaginary parts of the matrix elements pio(t) and poi(t), which corresponds 
to interference terms, decay more slowly then the element pn(t). 
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3.2. Small cavity 

For a finite cavity, the spectrum of eigenfrequencies is discrete, Aui is large, and so 
the approximation made in the case of large cavity does not apply; no analytical result 
can be obtained for foo{t), given by Eq. (I2"2l) . in this case. However, if the cavity is 
sufficiently small, the frequencies fl r can be determined as follows [8] . In terms of the 
small dimensionless parameter 

5 = — = — , 56 
Au TIC 



cotM) = - + -(l-^), (57) 



Eq. (113]) is rewritten as 

9 1 / _ ttu 2 5" 

where 9 = Q/Au. With 5 <C 1, which corresponds to R <C irc/g (a small cavity), as 
plots of both sides of Eq. (l5Tj) as functions of 9 are drawn, one sees that, for k — 1, 2, . . ., 
the solutions of Eq. ( |57l) can be written as 9k ~ & + with £)t< 1; then, expanding 
cot(7refc) for small e^, one finds 

n k *AJk + ^)= 9 - + (58) 



7r/c/ 5 \ Tlk J 

If we further impose that 5 < 2g 2 /nu 2 , a condition compatible with 5 <C 1, then f2 is 
found to be very close to lj, that is, 

O « Q (l - ^ ) . (59) 



To determine foo(t), we have to calculate the square of the matrix elements (£q) 2 
and (£°) 2 . From Eq. (fl~5l) . using that ojk = kAu = kg/5, we find (to first order in 5) 

(t° k ) 2 - ^ (^) 2 - (60) 
Now, considering the orthonormalization condition (111) , we determine (tp) 2 as 

(0 2 -(i+7E r j =( 1 + H"' (61) 

where we have used that £(2) = YlT=i ^~ 2 = 71 " 2 /6. We thus obtain, for sufficiently small 
cavities (8 <C 1), 

~' _, .n-iiw ! '' ( 

fe=l 

Notice that, if we had calculated (t[j) 2 by approximating directly Eq. <HM for small 5, 
instead of using the normalization condition (j3J), /oo(£) would not satisfy the requirement 
\foo(t)\ 2 < 1 with the equality occurring when t = 0. Also, it is worth mentioning that 
the approximation we made holds independently of the strength of the coupling between 
the particle and the environment and, thus, it applies for both weak and strong limits. 



/oo(t) 



\ / i — i 



(62) 
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Substituting foo(t), given by Eq. (16"2~|) . into Eqs. (J3HD (EHD leads to the elements of 
the reduced density matrix for the case of a small cavity. Distinctly from the case of 
a very large cavity, now, these matrix elements do not tend to zero as time evolves; in 
fact, they oscillate, never reaching zero. In particular, we find that 

k=l 



Pu{t) 



0(1 - 



6 [k+ TTk' 



t 



+ 



165 2 



7T 



k,l=l 



1 



COS 







t 



(63) 



Note that pn(t) is an oscillating, nonperiodic, function of time. A lower bound for p n 
can be found if one replaces the cosine terms in Eq. fl6"3"|) by —1; in fact, to order 5 2 , one 
finds 



Pu(t)>£ 



7T 



5 



2^2 
■ 7T 



(64) 



3 9 

This means that the particle state never reduces to the dressed-particle vacuum state; 
the particle system, maintaining itself in a mixed state, keeps exchanging quanta with 
the environment. 

The degree of impurity of this particle state at time t can be calculated using 
Eq. fl39|) or, equivalently, with 

D{t;0 = 2p n (t)[£-p n {t)]. (65) 

This quantity is plotted in Fig. [31 for the same values of u, g and £ as those used 
in Fig. [2] (for comparison), fixing a specific value of 5; one clearly sees the oscillatory, 
nonperiodic, behavior mentioned above. Therefore, in the case of a small cavity, for 
both weak- and strong-coupling regimes, the particle system is not dissipative; it keeps 
exchanging energy with the environment. For example, in the case of weak coupling, 
a physically interesting situation occurs when interactions of electromagnetic type are 
involved. In this case, we can take g = au>, where a is the fine-structure constant, 
a ~ 1/137. For a frequency Co ~ 2 x 10 11 /s (in the microwave band) and a cavity with 
R ~ 10~ 2 m, we have 5 ~ 0.016 and, so, the diagonal element pu(t) will never fall below 
the value ~ 0.87£. For strong coupling, g — (3u (j3 > 1), a similar result holds for a 
smaller cavity. The same kind of oscillatory behavior can be verified for the nondiagonal 
elements of the density matrix. In other words, the system confined in a small cavity is 
not dissipative, contrary to the case of a large cavity. 



4. Concluding Remarks 

We have presented in this paper a nonperturbative treatment of a quantum system 
consisting of a particle, in the harmonic approximation, coupled to an environment 
modeled by noninteracting oscillators. We have used renormalized coordinates in 
terms of which dressed states can be constructed. These states allow to naturally 
separate the system into the dressed (physically observed) particle and the dressed 
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D 
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. 1 

1 2 3 4 5 6 t 

Figure 3. Behavior of D{t\ £) as a function of t for 6 — 0.1, with Q = 1.0 and g = 0.5 
fixed (in arbitrary units), for some values of £: 0.3, 0.6 and 0.9 (dotted, dashed and 
full lines, respectively). 

environment by means of the conveniently chosen renormalized coordinates, q' and q'j, 
associated respectively to the dressed particle and to the dressed oscillators composing 
the environment. The dressed particle will contain automatically all the effects of the 
environment in it. Using this formalism we perform a nonperturbative study of the 
time evolution of a superposition of the ground state and the first excited state of the 
particle, the system being confined in a perfectly reflecting cavity of radius R. 

In the R —>■ oo limit, we find dissipation (damping) with dominance of the 
interference terms of the density matrix, in both weak- and strong-coupling regimes. 
For small values of R, the system is not dissipative, presenting stability. In both cases 
we have obtained results in agreement with expected behaviours. However, we do not 
simply recover well-known features with our method. For instance, as far as we know, 
the corrections to the exponential decay in Eqs. (|48p to ([51]) and from Eqs. (|52|) to 
(155]) have not been reported before in the literature. Also, the absence of dissipation 
in a small cavity can be understood as the analogue of spontaneous emission inhibition 
for an atom, which has a qualitative explanation when the lowest resonance cavity 
frequency is larger than the oscillator frequency. But also in this case we are able, with 
our method, to give rigorous expressions for the elements of the density matrix (in an 
analogous manner as has been done in [8] for the size of cavities ensuring stability of 
excited atoms). 

It is worth to stress that the renormalized coordinates and the dressed states, 
defined in Eqs. ( JTTjl and (jIBl) . allow an exact treatment of the problem, in which one 
completely avoids the use of perturbation theory. As already noticed in the text below 
Eq. (l20i) . the renormalized coordinates and dressed states are new collective objects, 
different from the normal coordinates Q and the eigenstates ([16]) . Distinctly from the 
eigenstates, our dressed states are all unstable, except for the ground dressed ({/^ = 0}) 
state. We assume that the dressed states (JTSl) are the physically meaningful states, 
instead of the ones written in terms of the bare coordinates g M . This can be seen as 
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analogous to the wave-function renormalization in quantum field theory, which justifies 
referring to q'^ as renormalized coordinates. Also, it is worth mentioning that the 
introduction of renormalized coordinates naturally ensures that the dressed vacuum 
state is stable, contrarily to the bare perturbative vacuum, since, by construction, it is 
identical to the ground state of the interacting Hamiltonian ([3]). Remember that the 
invariance of the ground state is due to our definition of renormalized coordinates given 
by Eq. (HZ}. 

As already explained in the Introduction, our formalism and results are restricted 
to a linear model, but it has the advantage of being nonperturbatively solvable. More 
realistic physical situations, however, would require a nonlinear coupling between the 
particle and the environment. It is nevertheless envisageable that such a study of 
complex nonlinear systems could be handled with our nonperturbative method. Indeed, 
the concept of renormalized coordinates has already been extended to a nonlinear model 
in Ref. [12j. Also, the generalization of the work presented in the present paper to the 
case of finite temperature, as well as the study of dressed coherent states, is in progress 
and will be presented elsewhere. 
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